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PACKET-ERROR  PROBABILITY  AN  ALYSIS  FOR  UNSLOTTED 
FH-CDMA  SYSTEMS  WITH  ERROR  CONTROL  CODING 


INTRODUCTION 


In  frequency  hopping  (FH)  systems,  a  code-division  multiple-access  (CDMA)  capability  can  be 
achieved  where  the  ra  patterns  take  on  the  role  of  codes.  Since  the  codes  usually  are  only  quasi- 
orthogonal  rather  than  truly  orthogonal,  frequency  hits  occur  and  lesult  in  loss  of  data.  In  the 
analysis  of  such  systems,  the  time-slotted  case  is  usually  considered  m  which  the  packet  length  is 
equal  to  the  siot  duration.  In  this  case  a  packet  transmission  is  sutgected  to  a  constant  number  of 
interfering  usera  throughout  its  duration.  When  Reed-Solomon  (RS)  error  control  coding  is  used  and 
all  frequency  hits  ate  assumed  to  result  in  symbol  errors,  the  packet-error  probability  can  be  evaluated 
straightforwardly  [1]  as 


Pr\packet  error  \  k  other  users] 


E 


1  -r+  1 


(1) 


where  is  the  symbol  error  probability  given  k  other  users  are  transmtting  simultaneously  over  the 
same  wideband  channel,  and  r  “  ~  *')/2j  is  the  error  correction  capability  of  the  RS(rt,v) 

code. 

In  this  report  we  consider  unslotted  systems  in  which  the  level  of  interference  varies  throughout 
the  transmission  of  the  packet  because  users  may  begin  their  transmission  at  any  time.  Consequently, 
the  symbol  error  probability  also  varies  throughout  the  packet  duratioii,  and  the  analysis  becomes 
much  more  difficult.  The  probkm  is  further  complicated  by  the'  dependence  of  symbol  errors  within 
the  packet,  since  the  interference  levels  experienced  by  the  symbols  of  a  packet  are  iiot  only  time 
varying,  but  also  they  are  dependent.  In  Ref  2.  Pursley  bounds  the  packet  error  probability  of 
unslotted  FH-CDMA  systems  in  which  the  interference  level  varies  over  the  packet  duration  by  the 
packet^rror  probability  of  the  system  with  the  maximum  interference  level.  Daigle  [3]  approximates 
the  packet  error  probability  of  an  unslotted  direct  sequence  CDMA  scheme  without  considering  the 
effect  of  error  control  coding. 

Often,  because  of  the  complexity  of  the  exact  analysis,  a  much  simpler  other-user  interference 
model  is  used  in  the  context  of  CDMA  systems.  In  this  model,  the  probability  of  packet  error  is  con¬ 
sidered  to  be  equal  to  one  if  the  number  of  other  transmissions  throughout  the  packet  duration  is 
greater  than  T,  a  threshold  le-'cl,  and  equal  to  zero  otherwise.  This  step  function  channel  model  has 
been  studied  in  Ref.  4  for  fixed-length  mesiages.  and  in  Ref.  5  for  exponentially  distributed  message 
lengths. 

In  this  report  we  present  an  exact  analysis  of  unslotted  FH-CDMA  systems  that  use  RS  coding. 
Fixed  length  packets  are  assumed.  The  computational  task  for  this  evaluation  is  enormous;  thus  it  has 
been  possible  to  evaluate  performance  only  for  small  codeword  sizes.  We  have  also  developed  upper 
bounds  and  close  approximations  to  the  packet-error  probability  that  permit  the  evaluation  of  more 

Manuicnp*  approved  January  5,  1989, 
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practical  systems.  We  have  demonstrated,  as  a  result  of  our  analysis,  that  the  bounds  that  are  based 
on  the  maximum  number  of  interferers  present  during  the  packet  duration  are  rather  loose,  and  that 
threshold  nrodels  do  not  provide  a  satisfactory  characterization  of  system  performance.  A  preliminary 
condensed  version  of  this  analysis  was  presented  in  Ref.  6,  and  a  more  detailed  development  can  be 
found  in  Ref.  7. 

SYSTEM  MODEL 

A  population  of  users  transmits  packets  containing  a  fixed  number  of  symbols  on  a  wideband  FH 
channel  consisting  of  q  orthogonal  narrowband  frequency  bins.'  One  ilf-ary  symbol,  representing 
IpgaAf  bits,  is  transmitted  per  hop.  The  hiding  patterns  are  assumed  to  be  generated  by  a  first-order 
Markov  process,  so  that  the  frequency  bin  for  each  hop  is  different  from  that  of  the  previous  hop,  but 
equally  likely  to  be  any  of  the  o^r  q  -  \  frequency  bins.  ' 

This  system  is  completely  asynchronous.  It  is  asynchronous  at  the  packet  level  in  the  sense  that 
packet  transmission  may  begin  at  any  time;  consequently,  packets  may  overlap  for  a  portion  of  their 
duration.  It  is  also  asynchronous  at  the  hop  level;  thus  fre^ency  hits  may  be  present  for  only  a  por¬ 
tion  of  the  duration  of  the  symbol.  However,  all  frequency  hits  are  assumed  to  result  in  symbol 
errors,  even  if  the  interference  is  present  for  only  a  small  fraction  of  the  symbol  duration.  The  result¬ 
ing  symbol  error  probability  given  that  k  other  users  are  simultaneously  transmitting  over  the  channel 
is,  as  discussed  in  Refs.  1  and  8, 

p*  =  1  -(1  -2/?)*(l  -po).  (2) 

where  po  is  the  probability  of  symbol  error  in  the  absence  of  other-user  interference,  i.e,,  the  symbol 
error  probability  caused  by  background  noise. 

Packet  sizes  are  chosen  so  that  each  packet  can  be  encoded  as  exactly  one  RS  codeword.  Each 
packet  consists  of  n  =  Af  A/-ary  symbols.  An  RS(n,  v)  code  can  correct  any  pattern  of  no  more  than 
T  =  [(n  -v)/2j  symbol  errors  in  a  codeword.*  We  have  considered  extended  RS  codes  of  rate 
1/2,  i.e.,  (n,n/2)  codes,  which  are  capable  of  correcting  any  pattern  of  no  more  than  r  = 
errors. 

PACKET-ERROR  PROBABILITY;  EXACT  ANALYSIS 

To  determine  the  exact  packet-error  probability  of  the  unslotted  FH-CDMA  system  with  RS  cod¬ 
ing,  it  is  necessary  to  determine  the  time-varying  level  of  other-user  interference  present  throughout  a 
packet’s  transmission.  We  consider  the  case  when  N,  the  total  number  of  other  active  channel  users 
during  the  transmission  of  a  particular  packet,  is  given.  To  describe  this  time-varying  interference 
process,  we  partition  the  interferers  into  two  groups:  “initial  interferers"  and  “final  interferers.” 
Consider  a  given,  user,  user  1 ,  whose  packet  transmission  begins  at  time  f o  ends  at  f o  +  T, 

where  T  is  the  fixed  packet  duration.  When  this  transmission  started,  li  packet  transmissions  were 
already  in  progress.  These  /,  channel  users  are  the  initial  interferers.  Each  of  these  will  terminate 
transmission  during  the  interval  (fo,to  T).  Now,  during  jser  I’s  packet  transmission  If  users 
begin  transmission  of  their  packets.  These  are  the  final  interferers,  who  will  still  be  transmitting  when 
user  1  ends  its  transmission  (see  Fig.  1).  Clearly,  /,  +  If  =  N.f 


•Undcfrctcd  codeword  error  probability  is  less  than  1/v?,  which  is  negligible  in  miny  applications. 

tit  is  assumed  that  (he  interfering  packets  cannot  arrive  iin«***«neousiy  with  the  tagged  packet.  Thus  an  interfering  packet  falls 
•in^mbigiKnish-  into  cifh-T  »hc  class  of  initial  or  bnal  intcrlcrers. 
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Now  define  the  departure  state  J  =<  (Ji,  j2,  ,  j^),  where  Ji  denotes  the  number  of  initial 

n 

interferers  that  terminate  transmission  during  the  /th  symbol  duration.  Note  that  ^  ji  •  U.  Like- 

_  t  “  1 

wise  define  the  arrival  state  k  ^  (k\,  ki,  . . . ,  k„),  where  ki  denotes  the  number  of  final  interferers 

that  begin  transmission  during  the  /th  symbol  duration.  Again,  we  have  ki  ^  If .  For  example, 

for  the  case  of  n  =4  symbols  per  packet,  /V  =  10  other  channel  users,  and  1/  =•  3  interferers  already 
transmitting  when  user  I’s  transmission  is  started,  one  possible  departure  state  is  (0, 1,0,2).  That  is, 
one  of  the  /,-  interferers  finishes  transmitting  during  the  second  symbol  of  user  I’s  packet  transmis¬ 
sion,  and  the  remaining  two  interferers  finish  during  the  last  symbol.  In  this  example.  If  is  equal  to  7, 
and  one  possible  arrival  state  is  (3,2,1, 1).  That  is  kj  =  3  interferers  start  transmitting  during  the 
first  symbol  of  user  I’s  transmission,  k2  =  2  start  during  the  second  symbol,  ^3  =  1  starts  during 
the  third  symbol,  and  k^  =  1  interferer  starts  during  the  last  symbol.  Note  that  numerous  possible 
departure  and  arrival  state  descriptions  represent  the  departure  of  the  the  If  initial  interferers  and  the 
arrival  of  the  If  final  interferers.  Given  If  and  N,  these  states  are  determined  by  an  exhaustive  search. 

We  assume  that  each  initial  interferer  is  equally  likely  to  terminate  transmission  during  any  sym¬ 
bol;  and  similarly,  each  fmal  interferer  is  equaUy  likely  to  start  transmission  during  any  symbol.  The 
probability  of  the  departure  state  j  and  the  probability  of  the  arrival  state  k,  given  If  and  N,  are  then 
given  by  the  multinomial  distributions. 


pr\j  I  If,  m  = 


Ii\{\/nf 
jl'Jl'-  .Jn'- 


Fr\k  I  If,  N]  = 


If\{\/n)‘’ 

kf\k2\...k„\' 


(3) 


Let  us  <)efine  the  interference  state  x  =  {x\,  xi,  . . .  ,  x„)  where  Xf  denotes  the  total  number  of 
interferers  present  during  the  /th  symbol  transmission.  The  value  of  Xf  is  determined  by 

i-\  I 

Xf  =  If  -  jy  +  ky  I  =  I,  2,  ....  n  (4) 

V  =•  1  V  «  1 

since  the  total  number  of  interferers  present  during  the  /th  symbol  is  equal  to  the  total  number  of  ini¬ 
tial  interferers  /,  less  the  number  of  initial  interferers  that  terminated  transmission  during  the  first 
(/  -  1)  symbols,  plus  the  number  of  final  interferers  that  arrived  during  the  first  /  symbols.  For  sim¬ 
plicity  we  consider  that  an  interferer  is  present  during  the  entire  symbol  in  which  it  starts  or  ter¬ 
minates  transmission,  even  though  we  do  not  require  symbol  synchronization.  This  is  consistent  with 
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our  earlier  assumption  that  all  hits  result  in  symbol  error  even  if  the  overlap  is  a  small  fraction  of  the 
symbol  duration.  Thus,  X[  is  actually  the  maximum  number  of  interferers  present  diuring  the  fth 
symbol^  In  the  above  example  the  interference  state  x  that  corresponds  to  the  departure  and  arrival 
states,  j  =  (0,1,0,2)  and  k  »  (3,2,1, 1),  is  J  »  (6,8,8,9).  Note  that  many  different  j  and  k  pairs 
may  combine  to  produce  the  same  interference  state  x. 

Now  the  probability  of  thejnterference  state  x  is  determined  from  the  probabilities  of  all  possi¬ 
ble  departure  and  arrival  states  j  and  k  that  combine  to  produce  x.  Given  the  state  x.  If,  and  N,  we 
have 

Pr[x\  It,  AfJ  =  2  PrO\  It,  N]  /V|itl  /„  N],  (5) 

litE. 

wherejE^  is_the  event  that  J  and  k  satisfy,  Eq.  (4).  The  probabilities  of  the  departure  and  arrival 
states  j  and  k,  given  /,■  and  N,  are  determined  from  Eq.  (3). 

Markov  Analysis 

A  Reed-Solomon  (n,v)  code  can  correct  any  pattern  of  i»  more  than  r  -  [(n  -  v)/2j  n-ary 
symbol  errors  in  a  codeword.  Thus  the  probability  of  packet  error  given  the  interference  state  x  is 
equal  to  the  probability  that  there  are  more  than  r  symbol  errors: 

n 

Pr\pkt  error  \  x)*'  =  '  Prlk  symbol  errors  \  x].  (6) 

*  ”  T+ 1 . 

A  Markov  analysis  approach  is  used  to  determine  the  probability  of  k  symbol  errors  given  x,  N,  and 
It.  We  consider  each  symbol  in  the  packet,  one  by  one,  starting  with  the  first  symbol.  As  each  sym¬ 
bol  is  considered,  we  determine  the  probability  that  the  symbol  is  in  error.  We  consider  the  probabil¬ 
ity  of  the  number  of  symbol  errors  in  the  packet  thus  far  as  the  number  of  symbols  considered  is 
increased  from  iv  to  w  +  1  for  1  <  w  <  n  -  1.  Thus  we  define 

+  (7) 


Pr{h  symbol  errors  in  the  first  w  +  1  symbols  |  I  symbol  errors  in  the  first  w  symbols]. 

Obviously,  the  only  possible  trwsitions  from  I  symbol  errors,  when  one  more  symbol  is  considered, 
are  to  /«  =  /  and  h  =  I  +  1.  Now  the  transition  from  /  to  /  +  1  symbol  errors  occurs  if  the  (w  +  l)th 
symbol  is  in  error.  Given  the  interference  state  x,  the  probability  of  the  (w  +  l)th  symbol  being  in 
error  is  given  as  in  Eq.  (2), 

Ph.  + 1  (/  + 1  I  /)  =  PrKw  +  l)th  symbol  in  error  |  x) 

=  I  -  (I  -  po)  (8) 


•Note  that  "pkt"  is  used  as  an  abbreviation  for  "packet”  througbout  this  report. 
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aod  thus, 

P*,+i(/|  /)  PrKw  +  l)th  jymho/  correct  f  x} 

•  d -2/gf-*'(l -Po).  <9) 

Now  define 

Pw(^)  "  symbol  errors  in  the  first  w  symbols  |  3r) .  (10) 

where  w  =  1,  2 . n  and  h  »  I,  2,. . n.  Thus, 

P^ih) PUh\  h)p^.i(h)  +  PUh\  h-l)p^.i(h-l)  (11) 

with  the  initial  condition  po(0)  =  1.  The  distribution  for  Pwih)  is  then  evaluated  recursively  until  we 
obtain 

p„(k)  =  Prlk  symbol  errors  \  x) ,  (12) 

from  which  it  follows  that  . 

Prlpkt  error  |  I)  =  Y,  />"(*)•  (13) 

*-r+l 

The  probability  of  packet  error  given  N  and  /;,  is  given  by 

Pr\pkt  error  |  //,  N)  =  '£iPr\pkt  error  \  xyPr\^\  /.,  N) 


=  E  E  /»"(*)  ^''(3?  I  //.  N\.  (14) 

I  . 

We  assume  that  /,  is  uniformly  distributed  between  0  and  N*.  Therefore, 

Pr\pkt  error  \  N\  =  E  E  E  ^1-  (13) 

-0;  L*-t+1  J 

The  computational  task  for  this  performance  evaluation  is  enormous.  Two  methods  were  used 
for  this  performance  evaluation.  The  first  method,  which  did  not  luc  Tiemory,  required  more  CPU 
time  than  the  second  method  that  used  memory.  Howev(  ant  of  memory  required  by  the 

second  method  for  packet  lengths  greater  than  five  symb<  :s  per  packet  was  prohibitive.  Calculations 
were  made,  using  memory,  for  packet  lengths  of  n  =  4  and  5,  and  without  memory  for  n  =  4,  5,  6, 
7,  8.  Table  1  summarizes  the  amount  of  CPU  time  required  for  these  computations.  Note  that  the 
packet  length  of  n  =  6  does  not  correspond  to  a  Reed-Solomon  code.  Calculation  of  the  performance 
for  this  packet  length  was  performed  with  an  arbitrary  value  of  r  to  establish  the  additional  amount  of 

•This  assumption  is  consistent  with  Poisson  arrival  suLstics. 
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CPU  time  required  when  the  packet  length  is  increased  by  one  symbol.  When  the  packet  length  is 
increased  by  one  symbol,  a  conservative  time  factor  was  detennined  from  the  computations  with  and 
without  memory.  The  CPU  time  required  if  infinite  memory  were  available  was  projected  based  on 
the  amount  of  time  required  for  n  4  and  5,  and  this  conservative  time  factor.  Note  that  even  if 
infinite  memory  is  assumed,  the  exact  perfortrumce  evaluation  of  the  common  packet  length  of 
n  »  16  symbols  per  packet  is  not  feasible.  Thus  efforts  were  made  to  develop  an  upper  bound  and 
approximation  to  the  packet-error  probability. 

Table  1  —■  CPU  Time  Required  for 
Packet  Lengths  of  »  »  4,  5,  6,  7,  8, 
and  16  Symbols  with  or  Without 
Memory.  Projected  CPU  Time 
Requir^  Is  Given  for  Computations 
with  Memory. 


CPU  Time  Required  (yr:h:  minis) 

A 

Memoiyless 

With  Memory 

~ 

00:00:00:59 

00:00:00:13 

5 

00:00:05:60 

00:00:00:58 

6 

00:00:24:40 

00:00:04:24* 

7 

00:01:52:24 

00:00:19:45* 

8 

00:05:46:35 

00:01:28:51* 

16 

25:00:00:00* 

‘Projected  CPU  time  required 


UPPER  BOUND  ON  THE  PACKET-ERROR  PROBABILITY 

All  upper  bound  on  the  packet-error  probability  was  developed  to  permit  the  evaluation  of  sys¬ 
tem  performance  when  practical  packet  sizes  are  used.  The  approach  of  the  upper  bound  is  described 
in  this  report.  The  transmission  interval  was  partitioned  into  halves,  quarters,  or  eighths,  each  con¬ 
taining  clusters  of  n  /2,  n  /4,  or  n  /8  symbols  respectively.  It  is  assumed  that  the  interference  level  is 
constant  over  the  duration  of  the  cluster  and  equal  to  the  maximum  level  of  interference  actually 
experienced  by  any  symbol  in  the  cluster.  That  is,  this  approximation  mixes  our  exact  analysis 
approach  with  the  approximation  used  in  Ref.  4.  Note  that  the  smallest  partition  interval  u  .cd  was 
eighths,  since  the  computation  of  the  exact  analysis  indicated  that  the  computational  limit  was  eight 
symbols  per  packet.  The  number  of  clusters  in  the  packet  for  a  given  partition  is  defined  as  b,  where 
b  =  2,  4,  or  8.  Each  cluster  then  consists  of  n/b  symbols.  The  new  interference  state 
X  =  (i|  ,^2.  •  •  describes  the  maximum  number  of  possible  interferers  in  the  new  clusters  of  the 
packet,  where  Jc,  denotes  the  maximum  number  of  interferers  in  the  ith  cluster  of  the  packet. 

We  must  recognize  that  this  new  interference  state  x  is  simply  tl  e  original  interference  state  x  of 
the  packet  with  b  longer  symbols.  The  length  of  the  new  symbols  is  equal  to  the  length  of  n  lb  old 
symbols.  With  this  knowledge,  we  realize  that  the  probability  of  x  with  n  symbols  per  packet,  and 
n  lb  symbols  per  cluster,  is  equal  to  the  probability  of  x  with  b  symbols  per  packet,  which  is  given  in 
Eq.  (5). 
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A  Markov  analysis,  similar  to  that  used  for  the  exact  analysis,  has  been  developed.  Here,  for  a 
given  partition  of  the  packet,  b  clusters  {b  »  2,  4,  or  8),  containing  n  lb  symbols  each,  we  consider 
each  cluster  of  the  packet,  one  by  one,  starting  with  the  first.  As  each  cluster  is  considered,  we 
determine  the  probability  that  I  symbol  errors  are  in  the  cluster.  We  consider  the  probability  of  the 
number  of  symbol  errors  in  the  packet  as  the  number  of  clusters  considered  is  increased  from  v  to 
V  +  1  for  1  <  V  <  n  -  1.  Thus  we  define 

/»»+i(A|0-  (16) 

Prlh  symbol  errors  in  the  first  v  +  !  clusters  |  i  symbol  errors  in  the  first  v  cbisters). 

The  only  possible  transitions  from  i  symbol  errors,  when  one  more  cluster  containing  n  lb  symbols  is 
considered,  are  to  h  =  i  +  /  symbol  errors,  where  I  =*  0,  I,...,  n lb.  Now  the  transition  fitjm i  to 
I  +  /  symbol  errors  occurs  if  there  are  /  symbol  errors  in  the  v  +  Ith  cluster.  Given  the  new 
interference  state  x,  this  is  given  by 

Pv+i(/  +  l  I  i)  =  Pr\l  symbol  errors  in  the  (v  +  l)th  cluster]  x] 

=  1^”/  Px.il  -p,,)"'*-',  (17) 

where  is  the  probability  of  symbol  error  given  x,  interferers  given  in  Eq.  (2).  As  before,  let  us 
define 

Pv(A)  =  Pr[h  symbol  errors  in  the  first  v  clusters  |  x) ,  (18) 

where  v  =  1,  2,  . . .  ,  b  and  h  ==  1,  2, . . .,  n.  Now, 

Pvib)  =  J^Pyih  \  h  -  k)py.iih  -  k)  (19) 

t=o 

with  the  initial  condition  po(0)  =  1,  The  distribution  for  Pvih)  is  then  evrJuated  recursively  until  we 
obtain 

Pb(k)  =  Pr{k  symbol  errors  in  b  clusters  |  x) 

=  Pr[k  symbol  errors  in  packet  |  xj ,  (20) 

from  which  we  can  form 

Pr\pkt  error]  xj  =  Pbik)-  (21) 

k—T’^l 
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The  probability  of  packet  error  given  It  and  N,  is  thus  bounded  by 

Pr\pkt  error  |  4,  S  J^Prlpkt  error  \  x]Pr^  |  Ii,  IV) 


1  J 


iV{x|  I„  N], 


and  thus 


Pr\pk$  error  \  N]  s  £  £  £  p»(*)  Pr^  |  I„  N] ,  03) 

since  4  >>  assumed  to  be  uniformly  distributed  between  0  and  N. 

APPROXIMATION  OF  PACKET-ERROR  PROBABILITY 

An  approximation  to  the  packet-error  probability  can  also  be  derived  to  ^rmit  the  evaluation  of 
the  packet  error  probability  for  longer  packet  lengths.  In  the  exact  analysis,  for  a  given  4  and  N,  we 
evaluated  the  probabilfty  of  packet  error  for  each  possible  inerference  state  x  for  the  transmitted 
packet.  Averaging  over  all  possible  x  produced  the  packet-error  probability  given  4  and  N.  For  the 
approximation,  given  4  and  N,  we  determine  the  expected  value  of  symbol-error  probability  for  each 
symbol  of  the  transmitted  packet.  Then,  irKorrectly  asssuming  that  the  symbol  errors  are  indepen¬ 
dent,  we  evaluate  the  packet-error  probability  when  the  symbol-error  probability  for  each  symbal  is 
given  by  its  expected  value. 

As  before,  we  partition  its  N  other  charuiel  users  into  fi  initial  interferers  and  If  fitud  in'xrfer- 
ers.  Starting  with  the  initial  interferers,  we  consider  each  symbol  in  the  packet  starting  with  tlie  first 
symbol.  Let  us  define 

ip  =  [n’tmber  of  initial  interferers  in  the  pth  symbol],  (24) 

The  probability  that  /,  initial  interferers  are  in  the  first  symbol  is  equal  to  1,  i.e.,  /’rjii  =  4)  =  l-O- 
As  in  the  exact  analysis ,  we  consider  an  interferer  to  be  present  during  the  entire  symbol  in  which  it 
starts  or  termiiiates  transmission;  thus  there  must  be  4  initial  interferers  in  the  first  symbol.  Now,  the 
probability  that  4  initial  interferers  are  in  the  second  symbol  is  the  probability  that  k  of  the  4  ini¬ 
tial  interferers  ended  their  transmission  during  the  first  symbol.  That  probability  can  be  expressed  as 

Fr(/2=4-k|i,  =4)=k  (25) 

rJ  L  "  J 

Now,  in  general,  the  probability  that  l—k  initial  interferers  are  in  the  (j  +  l)st  symbol,  given  that  I 
initial  interferers  are  in  the  yth  symbol,  is  the  probability  that  k  of  I  initial  mterferers  ended  their 
transmission  during  the  yth  symbol,  that  is. 


Prlij  .i  =  !  -k  \  ij  =  l]  = 


n  -  j  +  I 


-7  +  1 
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Now  the  probal)iUty  that  k  initia]  interfeten  are  m  the  (/  -f  I)th  symbol  gives  If,  and  AT  can  be  deter* 
mined 


i,-k 

Prlij^i  -  *  I  /,,  N)  -  2;  PrUj*i  “  *  I  Ij  -  k  +  m]  +  «i  |  /,,  AT)  (27) 

mmO 

aiki  the  initial  condition  PrUi  «  /J  ■■  1.0.  This  yields  the  distributkn  of  initial  interfeiers  in  each 
symbol  given  If,  and  N. 

To  obtain  the  distribution  of  final  interfeiers  for  each  synmrf,  we  again  consider  each  symbol 
in  the  packet  starting  with  the  first  symbol.  Let  us  define 

fp  *  {number  of  final  interferen  in  the  pth  symbol].  (28) 

Now,  the  probability  that  no  final  interfeiers  are  transmitting  prior  to  the  first  symbol  is  equal  to  1.0. 
We  call  this  the  probability  that  no  fuul  interferers  are  transmitting  in  the  0th  symbol,  i.e., 
Pr[fo  0|  =  1.0.  Now,  the  probability  that  k  Onal  interferers  are  in  die  first  symbol  is  the  proba* 
bility  that  k  of  the  If  final  interferers  started  transmission  during  the  first  symbol.  This  can  be 
expressed  as 


Prill 


k\  fo 


r  ’ 

k 

r 

n  - 1 

k 

n 

L  J 

L » J 

(29) 


In  general,  the  probability  that  I  +  k  final  interferers  are  in  the  (/  +  l)th  symbol,  given  that  /  final 
interferers  are  in  the  yth  symbol,  is  the  probability  that  it  of  the  //  -  k  final  interferers  started 
transmission  during  the  (j  +  l)th  symbol,  that  is. 


Prifj^i  =  I  +  k  \  fj  ^  I,  N,  It] 


The  probability  that  k  final  interferers  are  in  the  (j  -1-1)01  symbol,  given  It  and  N,  can  be  determined 
by 


+  ,  =*|  /„M  = 

i:  Prifj  ,1  =  k\fj  =  k  -  m,  S,  It]  Pr[fj  =  *-« |  If,  N]  (31) 

m  «  0 

and  the  initial  condition  Pr(/o  =  0)  =  1.0.  This  yields  the  distribution  of  final  interferers  in  each 
symbol  given  /,  and  N. 

The  distributions  of  initial  interferers  and  final  interferers  in  each  symbol  are  independent  (given 
Ij  and  A);  thus  the  distribution  of  the  total  number  of  interferers  in  each  symbol  can  be  obtained  by 

Pr{tj  =k  I  Ii,N]  =  £  Pr[ij  =  m  |  It.N]  Pr{fj=k-m  \  l„  N]  (32) 

rft  “  0 

where  tj  is  defined  as  the  total  number  of  interferers  in  the  yth  symbol. 
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We  can  now  determine  the  expected  >ahte  of  the  probability  of  symbol  error  of  the  ^tfa  symbol, 

Py  ■  2  Pr\xymbol  error  |  m  btterferen]  Pr[tj  -  m).  (33) 

m-O 

Let  us  defme  an  average  symbol  errcr  probability  state 

P^(puh . h),  (34) 

where  p,  is  the  expected  value  of  the  probability  of  symbol  error  of  the  /th  symbol.  Assuming 
independence  of  symbol  errors,  we  can  determine  the  probability  of  packet  error  given  the  average 
symbol  error  probability  state  p  as  in  the  exact  analysis.  As  mentioned  above,  this  assumption  is 
incorrect;  since  the  symbol  errors  are  not  independent,  this  results  in  only  an  approximation  to  the 
packet-error  probability.  Now,  we  need  to  consider  all  possible  occurrences  of  k  errors  in  the  code¬ 
word  and  form  the  error  state  e  “  (ei.ej, . . .,  e*),  and  the  correct  symbol  state 
d  «  . .  .,dn  ..h)  where  denotes  the  symbol  number  of  the  ith  incorrect  symbol,  and  d( 

denotes  the  symbol  number  of  the  ith  cotrect  symbol. 

Now,  the  probability  of  ?,  given  N,  /(,  and  p,  is  given  by 

Prl?  I  N.  A.  P!  -  n  K  n*(‘  -  (35) 

M  9  I  y  •  I 

and  thus  the  probability  of  packet  error  given  N  and  /,  can  be  approximated  by 
/•rip*/ errorl  )V./,1  -  £  , 

=  z:  E  -pj.)  .  (36) 

*-r+l  f,<r,  [“’■I 

where  is  the  set  of  all  ?  states  corresponding  to  k  symbol  errors,  and  the  RS  ctxle  used  has  an 
error  correction  capability  of  r  symbol  errors.  Now, 

Mpkt  error  I  N\  =  E  L  E  11  K  "O  ('  “  h)  .  (37) 

-  0  *  »  r*  I  -  I  V  -  1 

since  all  partitions  of  N  arc  equally  likely;  that  is.  /,  is  assumed  to  be  uniformly  distributed  between  0 
.ind  N  (again  consistently  with  the  assumption  of  Poisson  arrivals). 

Another  approximation  to  the  packet-error  probability  can  be  developed  by  calculating  the 
expected  number  of  interferers  transmitting  during  each  symbol  from  the  probability  distribution  of 
t!ie  total  number  of  interferers  in  each  symbol.  That  is.  given  /V|f,  -  k  |  /..At.  calculate. 

,v 

Xj  =  E  ">  ^'•1'/  ^  "*1  A.  'V|.  (.18) 

m  ■»  () 
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Then  the  probability  of  symbol  error  of  the  Jth  symbol,  given  die  expected  mimber  ci  interferers  ia 
the  symbol,  could  be  determined  as, 

Pj  ■  Pr{symbot  error  \  xjinterferers]  (39) 


-  I  -  0  -po). 

Substituting  pj  for  pj  in  Eq.  (37)  yields  the  second  approximation. 

SYSTEM  PERFORMANCE 

Packet-error  probability,  given  N,  was  evaluated  for  several  rate  1/2  RS  codes  as  N  was  varied 
between  I  and  10.  In  ail  cases  the  number  of  frequency  bins  wu  f  SO,  and  a  noiseless  channel 
(i.e.,  po  **  0)  was  assumed.  Exact  performance  results  were  obudned  only  for  the  RS(4,2)  and 
RS(8,4)  codes.  Upper  bounds  and  approximations  were  obtained  for  these  codes  as  well  as  for 
RS(  16,8)  and  RS(32. 16)  codes. 

The  RS(4,2)  aitd  RS(8,4)  codes  have  error-correction  capabilities  of  r  >■  I  and  2  symbol  errors 
respectively.  These  are  not  realistic  packet  lengths  since  they  permit  the  transmission  of  only  very  lit¬ 
tle  information  and  since  the  upper  bound  on  undetected  codeword  error  probabilities  for  the  RS(4,2) 
and  RS(8,4)  codes  are  equal  to  1.0  and  O.S  respectively.  However,  the  computation  of  the  exact 
packet-error  probability  of  the  system  with  the  short  packet  lengths  it  very  usefol;  it  allows  com¬ 
parison  with  the  upper  bound  and  approximation  to  the  packet-error  probability.  Computational  limits 
prevent  the  exact  evaluation  of  the  packet-error  probability  for  longer  packet  lengths  at  (!iscus.sed  ear¬ 
lier.  The  upper  bound,  with  n  /h  symbols  per  cluster.  b»  2.  4,  and  8,  was  computed  for  the  ,RS 
codes  mentioned  above,  as  well  at  the  RS(I6,8)  and  RS(32,16)  codes  having  error-correctiou  capabili¬ 
ties  of  r  *  4  and  ,8  symbol  errors  respectively. 

Figure  2  shows  the  results  for  the  unslotted  FH-CDMA  system  with  RS(4.2)  code.  Note  that 
the  approximation  is  er.trcmely  close  to  the  exact  packet-error  probability.  When  the  number  of 
interferers  transmitting  during  the  packet  N  is  greater  than  or  equal  to  3,  the  approximation  is  within 
5%  of  the  exact  packet-error  probability,  whih  it  is  within  1%  for  A/  a  7. 

The  upper  curve  in  each  of  the  figures  represents  a  slotted  system  in  which  all  N  interferers  are 
present  during  the  entire  packet  transmission. 

Figure  3  illustrates  the  performance  of  the  unslottcd  FH-CDMA  system  with  the  RS(8,4)  ctxle. 
As  before,  the  approximation  closely  resembles  the  exact  packet-error  probability.  When  A/  a  4  the 
approximation  is  within  10%  of  the  exact  packet-error  probability,  while  it  is  within  3%  when 

N  >  7. 

For  the  system  parameters  considered,  the  appniximalion  pttxvidcs  excellent  agreetnent  with  the 
cx:i(  f  resiilis  It  appears  that  the  proportional  difference  between  the  exact  packet-error  probability 
and  !tie  approximation  decreases  as  N  increases.  Also  note  that  the  packet-emir  probability  as  a  func¬ 
tion  of  N  docs  not  resemble  a  step- function,  thus  indicating  'hat  a  threshold  model  does  not  provide  a 
good  indication  of  system  performance,  as  demonstrated  in  Ref.  9. 

Figures  4  and  5  illustrate  the  upper  bound  and  approximation  of  the  performance  for  the  cases  in 
which  the  exact  performance  is  not  computable.  Figure  4  shows  the  performance  of  the  unslottcd 
FH  ('I)MA  system  with  the  RS(lb.8)  cixle.  Figure  .3  illustrates  the  approximation  and  upper  bound 
(if  the  packet-error  probability  of  the  system  with  the  RS(32.1ft)  code. 
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Pi«.  2  —  PiclHl-errer  probabUKy  wWi  RS<4,2)  codiag 
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It  is  of  interest  to  compare  the  packet-error  probability  of  the  system  with  time-varying  interfer¬ 
ence  to  that  of  the  system  with  constant  interference.  Table  2  contains  the  number  of  constant  interfer- 
ers  that  produce  the  same  probability  of  packet  error  as  10  interferers  whose  transmissions  start  or 
stop  within  the  packet  duration,  where  we  ignore  the  fact  that  this  number  must  be  an  integer.  The 
comparison  is  given  for  RS(4,2),  RS(8,4),  RS(16,8),  and  RS(32,16)  codes  and  9  -  50  and  100  fre¬ 
quency  bins.  For  packet  lengths  of  n  »  16  and  32  symbols,  the  approximate  packet-error  probability 
for  =  10  partial  interferers  was  used  for  the  comparison.  Note  that  this  number  appears  to 
approach  5  as  the  code  length  increases. 


Table  2  —  ConqNtrison  of  Time-Varying 
Interference  to  Cor^tant  Interference 


Constant  Interference  Level 

Corresponding  to  10  Partial  Interferers 

Code 

q  *  50 

q  -  100 

RS(4.2) 

6.1 

6.3 

RS(8,4) 

5.5 

5.6 

RS(16,8)* 

.  5.15 

5.15 

RS(32,16)* 

5.1 

5.1 

*Approxinuite  p*cket-«rror  probtbility  for  N  •  10 
partial  interferen  uitd  for  comparison. 


CONCLUDING  REMARKS 

The  packet-error  probability  performance  of  unslotied  FH-CDMA  with  RS  error  control  coding 
is  evaluated  exactly  for  small  packet  lengths.  Computational  limits  prevent  the  exact  evaluation  of 
packet-error  probability  for  longer  packet  lengths,  aivd  thus  an  upper  bound  and  approximation  to  the 
packet-error  probability  are  derived  and  computed.  The  upper  bound  calculations  are  also  limited 
(though  not  quite  as  severely)  by  computational  considerations,  but  the  approximation  is  easily 
evaluated  for  large  codeword  sizes. 

In  all  cases  considered,  the  packet-error  probability  as  a  function  of  the  number  of  other 
transmitting  users  Joes  riot  resemble  a  step  function.  Thus  detailed  models  that  reflect  channel  charac¬ 
teristics  and  coding  properties  are  needed  to  provide  an  accurate  evaluation  of  system  performance. 
Our  previous  observation  (9)  that  the  step-function  channel  model  docs  not  reliably  predict  perfor¬ 
mance  of  spread  spectrum  multiple  access  systems  is  further  confirmed  by  this  study. 
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